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Abstract 

We derive the mixed as well as unmixed lattice heavy meson chiral Lagrangian up to order 0(a 2 ), 
with Wilson and Ginsparg-Wilson fermions. We consider two flavor partially quenched QCD and 
calculate vector meson mass corrections up to order 0(a 2 ), including the corrections associated 
with the violation of rotational 0(4) symmetry down to the hypercubic group. Our calculations 
also include the one-loop, phenomenological contribution from the p — > nit decay channel. The 
final result is a chiral-continuum extrapolation formula with model dependent coefficients from 
which one can recover the physical p meson mass from the large amount of current lattice data. 
As a verification of our result, the chiral-continuum extrapolation formula is compared with that 
used in numerical simulations. 
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INTRODUCTION 



In the present day lattice simulations the quark masses used are much heavier than 
those in nature because of the increased computational cost at low mass. Instead of 
trying to reach to lighter quark masses one runs lattice simulations with heavier quarks 
and then performs a chiral extrapolation to the physical regime using the predictions of 
chiral perturbation theory (ChPT) 0, H|. However, it is not totally consistent to apply 
ChPT for extrapolation of lattice results because it is a continuum theory which doesn't 
contain information about the lattice spacing. That is why it is first required to extrapolate 
lattice data to the continuum limit and then perform the chiral extrapolation to physical 
regime. In order to consistently include lattice discretization effects, Symanzik's effective 
theory needs to be appli ed p , 0] • This approach is used, for example, in Refs. 0, 0, 0| and 
generalized in Refs. [H, H, llClj] to study the mixed lattice theories with Wilson fermions as sea 
and Ginsparg- Wilson fermions as valence quarks. Being relatively expensive to simulate, 
Ginsparg- Wilson fermions on the other hand, have an exact chiral symmetry even on 
the lattice. The latter allows one to exploit the good chiral properties of Ginsparg- Wilson 
fermions by considering them as valence quarks 1 and considering Wilson fermions as sea 
quarks 2 ■ The Lagrangian approach used here is necessary because of the large lattice 
spacing effects in unquenched simulations (which means that the 0(a 2 ) effects need to be 
accounted). It represents an efficient method by which to improve the numerical accuracy 
of the physical results. 

In this paper we calculate the vector meson mass corrections up to order 0(a 2 ) in par- 
tially quenched chiral perturbation theory (PQChPT) with graded symmetry SU(4\2), using 
Wilson and Ginsparg- Wilson fermions and the following power counting scheme: 

p 2 ~ a ~ m q (1) 

As the small parameters of our double expansion we will choose: e 2 = (aA,m q /A), so that 
m q < A C 1/a, (where A is a typical QCD scale). 

To calculate the corrections coming from lattice discretization one must formulate the 
corresponding lattice theory and match the new operators with those in chiral perturbation 
theory (ChPT). The standard procedure consists of the following parts: 



1. writing the Lattice QCD action for the fermions one is going to use; 

2. finding the corresponding Symanzik action up to a given order in a based on symmetry 
constraints from underlying lattice theory; 

3. performing a spurion analysis with further projection to ChPT. 

If we want to find the mass corrections to the vector mesons, we need to formulate ChPT 
for both pseudoscalars and vector mesons. ChPT for pseudoscalars with Wilson quarks was 
formulated, for example, in Ref.j?}. It was combined with Ginsparg- Wilson fermions as a 
mixed theory in Refs.jH, Continuum heavy meson ChPT was formulated in Ref.[l3| 



such valence fermions will have masses smaller than the valence quark masses accessible using Wilson 
fermions 
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and further developed for QChPT and PQChPT in Refs.|14| and [15|. An example of two- 
flavor PQChPT can be found in Ref . [l8l| . The main disadvantage of similar heavy vector 
meson ChPT formulations was that they didn't include terms which explicitly break 0(4) 
rotational symmetry down to the hypercubic group. Recent studies, particularly by Tiburzi 
in Ref. [loj . with Wilson and Ginsparg- Wilson quarks, showed explicitly the terms in the 
heavy baryon Lagrangian that break 0(4) symmetry, up to order 0(a 2 ). Following similar 
lines, we formulate heavy meson partially quenched ChPT (HM PQChPT) for Wilson and 
Ginsparg- Wilson fermions up to order 0(a 2 ), including chiral as well as continuum symmetry 
breaking. 

The paper is organized in the following way: in Sec. II we review the well developed 
PQChPT SU(A\2) Lagrangian for the pseudoscalar meson sector with mixed and unmixed 
Wilson and Ginsparg- Wilson fermions up to order 0(a 2 ); in Sec. Ill the extended heavy 
meson PQChPT Lagrangian is derived up to order 0(a 2 ), including the terms explicitly 
breaking 0(4) rotational symmetry down to the hypercubic group (this is actually the ad- 
justment of Ref. [10( applied to vector mesons); in Sec. IV we calculate the p meson mass 
corrections including phenomenological one-loop correction from p — ► tttt and in Sec.V the 
chiral-continuum extrapolation formula is derived and compared with the phenomenological 



expression in Ref. [21J, in terms of which the p meson mass was successfully extracted. 



PSEUDOSCALAR MESON SECTOR 

Let us denote the quarks of S77(4|2) as 

Q = (x,y,u,d,x,y) T , (2) 

where x and y label the valence fermions whose ghost partners (x, y) cancel their effects in 
the sea loops. Thus, the only sea loop contributions come from unquenched u and d quarks. 
Now, assume that m x = m y and m u = which is natural if we want to calculate the 
correction to the mass of the p meson. The mass and Wilson matrices l]| are then: 



m Q = diag(m x , m x , m u , m u , m x , m x ), (3) 
w Q = diag(w v , w v , w s , w s , w v , w v ) (4) 

where w v refers to valence and w s to sea quarks. Note that if quark is a Wilson fermion 
then: {wq)i = 1, while if it is GW then: (wq), = 0. This could be useful for the construction 
of mixed types of ChPT, as in Refs. |H 

The leading order (LO) Lagrangian, L 2 ~ 0(p 2 ,a,m q ), which includes chiral symmetry 
breaking from quark mass and lattice discretization, is given by the expression ( [8l [To| ) : 

f 2 

U = y < d^d^ > - \ m < m Q T) + mjjS > - a\ a < w Q ^ + w^E > (5) 

where the angular brackets stand for the super-trace over the flavor indices, / is the pion 
decay constant, A m ( a ) are low energy constants and £ = exp(2iQ/f) contains the matrix 
of meson fields, $ G S77(4|2). The next to LO (NLO) Lagrangian would be of order 
0{e A ) ~ 0(p 4 ,p 2 a,p 2 m q ,a 2 ,m 2 am q ), which is rather lengthy (see Ref. 0]). However, in 
calculating the vector meson mass correction the NLO pseudoscalar meson Lagrangian will 
give corrections of 0(e 5 ) in which we are not interested. We also explicitly omitted from 
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the Lagrangian the contribution from the rj (the SU(2) analog of it) which, because of the 
U(1)a anomaly, is integrated out from the theory as a heavy degree of freedom. However, 
double hairpin vertices involving the rj' were used to find disconnected propagators. 
The following results from Ref. are needed for our calculations: 

4 

m QQ' = J2 [ X mirn Q + m Q /) + a\ a {w Q + w Q >)} , (6) 

r _ i (g 2 + m 2 8ea ) _ i (r 

where raqqi is the mass of the pseudoscalar meson with quarks Q and Q' and m sea , m va[ 
and m m i X are mesons which consist of sea-sea, valence-valence and valence-sea quarks cor- 
respondingly. In our case: 

g 

m lea = -p l X mm u + d\ a } , (8) 
g 

m lal = p [Kn,m x + a\ a ] , 

Finally, G xy is the disconnected meson propagator, which includes all diagrams with double 
hairpin insertions and G vs is the connected propagator for the meson lines which consist of 
one sea and one valence quark. 



VECTOR MESON SECTOR 

Standard chiral perturbation theory for vector mesons has been formulated by Jenkins, 
Manohar and Wise in Ref. the quenched counterpart by Booth, Chiladze and Falk in 
Ref. |l4| and the partially quenched counterpart by C. K. Chow and S. J. Rey in Ref. 1511 . 
The study of vector meson masses within the context of ChPT 2 can be found in Ref. 17( • 
Here we will construct the partially quenched theory using mixed type of fermions and 
including the lattice discretization effects up to order 0(a 2 ). 

The heavy vector meson Lagrangian, up to order 0(e 4 ), consists of the following parts: 

L« ~ C(p), (9) 
~0(o,m 9 ), 
~ 0{a 2 ,m 2 q ,am q ). 

The Lagrangian is omitted here because it contains derivatives which can be eliminated 
using the equations of motion [20]. In the partially quenched (PQ) case the vector meson 



performing an expansion in terms of the momenta, quark masses and 1 /N c 



4 



multiplet (in the large N c limit) is described by a 6 x 6 matrix field (similar to the Goldstone 
meson sector): 

^=(*t*) ( 10 ) 

where V is a (4 x 4) matrix of qq states, V is (2 x 2) matrix which consists of qq states, 
\1/ and W are rectangular matrices involving fermionic bound states with qq or qq content. 
Under the G = iS77(4|2)l x SU(A\2) r x U{l)y graded chiral symmetry group, the heavy 
meson field transforms as: 

N^UN^Ul UeG (11) 

and under charge conjugation: 

CN.C- 1 = -Nf. (12) 

As in conventional ChPT (as well as QChPT [14] and PQChPT 15]), the vector meson 
multiplet will be treated as heavy, static source. We will also assume that the static vector 
meson propagates with a fixed four- velocity t> M , v 2 = 1 and interacts with soft Goldstone 
multiplets. The three polarization states of the vector mesons are perpendicular to the 
propagation direction (i.e. v ■ N — 0). 

Taking all of these into account, let us now proceed to construct the relevant Lagrangians. 
The LO heavy-meson chiral Lagrangian will consist of terms of order 0(p) ~ 0(e) which 
are: 

L kin = -%< Nl(v ■ D)N fl >-iA l <Nl>{v-D)<N ll > ~ 0(e), (13) 

L int = igx < Nl >< N U A X > v^ vX ° + h.c. + ig 2 < {JVj, N U }A X > v a e^+ (14) 
ig 3 < Nl >< N v >< A\ > v a e^ + i 9i < >< A x > v a e^ Xa ~ 0(e), 

where the following notations are implied: 

D,=d, + \v^ \, f^^'+w, = ~ (e^e f - e^e) (15) 

£ = M=i((tm (t + (m Q (). (16) 

Here we treat A\ as a small quantity, so the effects of this term can be reabsorbed in the 
other parameters by the field redefinition: 

N'» = N„ + y <N »>- 

This redefined field still transforms under the chiral group in the same way as N^. As we 
may notice, there are four-types of chiral invariant interactions in Eq. ()14j) (consistent with 
graded chiral symmetry) between vector meson multiplets and pseudoscalar mesons. In the 
large N c limit the coupling constants scale with N c as follows: g% ~ 4-, g 2 ~ 1, #3 ~ vw, 
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(74 ~ jr- However, according to Ref. the corresponding quark (and color) flow diagrams 
are not planar, so they need to be multiplied by an additional factor of 1/N C . 

Because we are interested in finding the one-loop correction to the flavor charged vector 
meson mass in isospin limit there will be no contribution from the terms proportional 
to g% and (73. These terms are responsible for quark flow diagrams with disconnected 
vector meson line, and because we require to have two different flavors at valence lines we 
cannot have this type of hairpins (as we don't consider any flavor changing interactions). 
The diagrams from the term proportional to g± are one hairpin like diagrams (due to 
coupling of the vector meson with the two flavor analog of r/). Sum of these diagrams 
with vacuum bubble insertions can be shown to vanish after the decoupling of flavor singlet 3 . 

The next to LO (NLO) Lagrangian - L 2 ~ 0(e 2 ), according to our counting scheme, will 
consist of terms of order 0{m q ) and 0(a). The first Lagrangian L mass ~ 0(m q ) is: 

— fi < N^Np > + m< Nj t ><N fl > + 

+ 01 < iVjiVp >< M > + a 2 < JVj >< >< M > 

+ X 1 (<Nl>< N^M > +h.c.) + A 2 < {JVj, N^M >~ G(e 2 ), (18) 

where the first term in Eq. p8j1 corresponds to residual mass p, of vector meson multiplets 
which is always possible to remove by a suitable reparametrization transformation (l6| . 
(through the choice of the reference momentum). Here we will use this freedom and assign 
Jx = 0. The second term is small because of Zweig's rule, therefore we treat it as the same 
order as the quark masses, 0(e 2 ) jlTf- The chirally invariant terms in Eq. ([18)1 proportional 
to a\ t 2 are contributing directly to the vector meson mass. For further simplicity and in order 
to match the conventional ChPT in Ref. |l3[ in continuum and large N c limits we will assign 
a\ t 2 = 0. The last two terms in L mass , proportional to Ai,2, correspond to SU(2) isospin 
breaking because of the quark masses. To match with the mass Lagrangian in Ref. [13| 
(Eq.24) we will take fi\ = Ai = throughout this paper reminding that this particular 
choice will be justified only by the final results. 

From the analysis above it follows that the flavor non-diagonal heavy vector meson 
propagator is: 

(v^n (19) 

v ■ k + xe 

where k^ = py — (My — /2)f M is its residual momentum, p v and My are the four-momentum 
and mass of the vector meson correspondingly. 

The second NLO term L a ~ 0(a) could be written in the following form: 

L a = ai (< M >< N^W + > + h.c.) + a 2 < {M, N fl }W + > (20) 



+ a 3 <Nl >< N„ >< W+ > + a 4 < N^N^ >< W+ > ~ C(e 2 ) 



3 For example, in Ref. 15] it is assumed that (71,3,4 = and gi = 0.75. In principle, our considerations are 
more general and we could leave these couplings as free parameters which are to be determined by further 
lattice simulations. 
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with the definition 



nA 2 

W + = — (fV + ^) > (21) 



where «i, 2,3,4 are undetermined low energy constants. Notice that the terms in L a break 
chiral symmetry in the same way as masses and that the dominant (with respect to N c ) 
term is one proportional to a 2 - The terms proportional to 0^3 don't give any contribution 
to the vector meson mass because (similar to gi t s) these interactions give hairpin insertions 
in vector meson line. 

Now, let us write the next to NLO (NNLO) Lagrangian of order 0{e i ). As we mention in 
the beginning the NNLO Lagrangian will consists of the following parts: L am ~ 0(am q ), 
L m 2 ~ 0(m? g ) and L a 2 ~ 0{a 2 ). 

Using the properties of M and W+ under chiral transformations the first term is: 

L amq =(3i<Nl >< A M X W+ X M > + fa < N^Nfj, >< W+XM> (22) 
+ (3 3 < JVj >< A M >< W + M > + p 4 (< >< N IM xM> +h.c) 

+ p 5 (< A+M >< N„ >< W+ > +h.c) + fa < N^Nfj, >< W+M > 
+ p 7 (< NlW+ X N^M > +h.c) + P 8 < {iVj, N^}W+ X M > 
+ (3 9 < {JVj, A M }M ><W+> + p 10 (< N} >< N„{W+, M} > +h.c) 
+ p 11 <{Nl,{N li ,W+}}M> 

Similarly, the second term has the structure: 

L m 2 = n < iVj >< JV„ >< M X M > + r 2 < Nl >< A M >< MM > (23) 
+ r 3 < N^N,j, X M >< M > + r 4 (< N^W+ >< N ll ><M> +h.c.) 
+ r 5 < >< MM > + r 6 < A^M >< N^M > 

+ r 7 < {Nl N„}M X M > +r 8 (< JVj >< {A M , M}M > +h.c.) 
+ r 9 <{Nl,{N„M}}M>, 

where Pi,.. .11 and r± i ... j g ~ 1/Ai are undetermined constants (A x is some energy scale). The 
terms r x 5 8 are 1/A C suppressed, according to Ref. [l7j . These Lagrangians give only tree 
level contributions to the mass at 0(e 4 ) which are proportional to am, and mj? correspond- 
ingly. 

The NNLO Lagrangian of order 0(a 2 ) will consist of four parts, the first part has the 
form: 

L% w = q 1< Nj 1 X Nn X W+ X W+ > + q 2 < A+ >< A M >< W+W+ > (24) 
+ q 3 < X W+ X W+ > + q A (< ^W+ X A M X W+ > +h.c.) 

+ q 5 < N^Nfj, X W+W+ > + q e < N}W+ >< N„W+ > 
+ q 7 < {At, N^W+ X W+ > +q 8 (< JVj >< {A^, W + }W + > +h.c.) 
+ q 9 <{Nl,{N„W + }}W + >, 

where gi,...,g ~ — 1/A. 
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In order to construct the second and third parts of 0(a 2 ) Lagrangian, taking into account 
the relation wq = 1 — wq (from fib]]), let us define: 

p + = \ [e f (^Q - vo Q )t + i{ WQ - w Q )e] (25) 

which transforms under G as UP + U^ . Using this spurion field the second part of the 0(a 2 ) 
Lagrangian will be: 

= 7o < NlN, > + 7l < iVj >< iV M >< P + > + 72 < >< P + > (26) 

+ 73 < {Nl N p }P + > +74 (< Nl >< N,P + > +h.c.) , 

where 7 o,...,4 ~ — a 2 A 3 /A\. Similar to Eq. (j23), the third part of 0(a 2 ) Lagrangian will be 
written in the form: 

If/ = 5t < JVj >< N p >< P + ><P + > + 5 2 < iVj >< N p >< P + P + > (27) 
+ 5 3 < N^ >< P + ><P + > + 8 A (< NlP + >< N p >< P + > +h.c.) 
+ 6 5 < N^ >< P + P + > + 6 6 < NlP + >< N P P + > 
+ 6 7 < {Nl N p }P + ><P + > + 5 8 (<Nl >< {N„ P + }P + > +h.c.) 
+ 5 9 <{Nl,{N„P + }}P + >, 

where <5i,...,9 ~ — a 2 A 3 /Ai. 

Finally, following Ref. [HJ], we write the last part of 0(a 2 ) Lagrangian, which explicitly 
break 0(4) symmetry, in the form: 

L c a2 = hvl < NlN v > + k 2 vl < Nl{N u , P + } > + k 3 vl < NtN v >< P + > 

+ k 4 vl < NlN, > + hvl < Nl{N„ P + } > + k % vl < NlN, ><P + >, (28) 

where fcj,...^ ~ — a 2 A 3 /Ax. Later we will see that for the vector meson at rest the mass 
corrections from lattice discretization up to 0(e A ) are of the form: M p (a) = M p (0) — a 2 AD, 
where D is some constant. 

VECTOR MESON MASSES UP TO ORDER 0{a 2 ) 

This developed heavy meson (HM) PQChPT with S77(4|2) graded symmetry allows us 
to calculate the corresponding vector meson mass correction. The vector meson mass in 
both chiral and continuum expansions can be written in the form: 

M 2 = [MM + + a^i) + afljji) + ...) 2 + + a^i) (29) 

where Mo(fi) is the renormalized vector meson mass in continuum and chiral limits, 
^rleiioop)^) * s ^h e ^ Tee oo P) contribution to the self energy of order 0(e n ) and jj, is the 
renormalization scale. At order 0(e 2 ) we have the following tree level contribution to the 
mass: 

a trle = ~^i m x — 2aA 2 [a 2 w v + a^Ws] . (30) 
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We will not have any contribution of order 0(e 3 ) simply because we don't have Lagrangian 
of this order, so cr\ ree ([i) = 0. 

Now, let us write the 0(e 4 ) correction as: 

= [«A 2 (Am x + Bm u ) + Am 2 x + Bml] 
Ai 

- a 2 A 3 (C + DvlNlK + DvlNlN,) , (31) 

where A, B, A, B, C, D and D are undetermined coefficients 4 . 

The contribution from the tadpole-like loops, with double hairpin insertion is: 



(4) 4a 2 aA 



A? 



w s + w v )F(m mix , //) - 2w v F(m va i, fi) - 2w v Am 2 ^ >F ^ riv ^ 1 ? W 



(32) 



2 

F(m, /i) = m 2 Zn— (33) 

where Am 2 = m 2 aZ — m 2 ea = -p(m x — m u ). The corrections aj^ p were calculated in Ref. (2*l| . 
including the loop contribution from p — > 7r7r (which is included phenomenologically) . Rather 
than repeating the expressions used in that analysis of data over a wide range of meson mass, 
which were based on finite range regularization, we show just the leading or next-to-leading 
non-analytic behavior. 



CHIRAL EXPANSION FORMULA 

The complete chiral expansion for the vector meson mass at rest is then: 
My = (M (/i) - Xi m x - aA 2 X2m u - Xi m l ~ Xz m u ~ fl A 2 [a 2 w v + a 4 w s ] + c b a 2 ) 2 (34) 



+ 



4«2aA 



2 



(w s + w v )F(m mix , fi) - 2w v F(m vah fi) - 2w v Am 



dF(m val ,[i) 
9m 2 val 



+ S tot (/i) 



where xi = 2A 2 + AaA 2 /Ai, X2 = B/A u xi = A/A 1 , X2 = B/A x and c b = —A 3 D' are 
new parameters as functions of the previous ones. The term in square brackets comes 
from tadpole diagrams and is the total one-loop correction from Ref. 21] which includes 
the 0fc2p(/i) contribution. If one fixes the masses of the sea quarks and introduces new 
dimensionful (reparametrized) extrapolation coefficients, then we will have: 

M v — ( c o ~ A 2 (a 2 w v + a A w s + X2 r m u )a + c b a 2 + (x + ca)m 2 + x m t) 2 (35) 

- ^JT~ {K + w v )F{m mix ) - 2w v F(m val ) - 2^Am 2 ^^ m ' ) | + S to '(/i), 

where m 2 = and Co, c, Xi an d X are new independent parameters which are functions 

of the old ones. 

Now, consider two different cases following from the theory: 



4 for example, A = (2/3 5 + 2f3 9 ) w s + (2/3 7 + A/3 W + Aj3 n ) w v and B = (A + + As + As) w s + 
(2/3 4 + 2/3 8 ) w v 
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1. when the sea and valence quarks are Wilson- type fermions then the expansion takes 
the form: 

My = (co - A 2 (a 2 + «4 + X2"m u )a + c b a 2 + (x + ca)ml + x m Vf (36) 

8a 2 aA 2 ( 2 aF(m^)\ tot 

F{m mix ) - F(m val ) - Am } + £ (//), 



2. when the valence quarks are Ginsparg- Wilson fermions and the sea quarks are Wilson 
fermions we have: 

M v = (cq - A 2 (a 4 + X2m u )a + c 6 a 2 + (x + ca)m 2 + x m Vf (37) 



4a 2 «A 



2 



■F(m mix ) + E tot (/i) 



The analysis in Ref.[21], which corresponds to the first case, was not sensitive to the existence 
of tadpole-like diagrams which are proportional to «2, so there is a good reason to believe 
that «2 is consistent with zero. Similarly, c and + Xi m u or «4 and X2 are also consistent 
with zero 5 according to the analysis in the same reference. Comparison with the numerical 
simulations in Ref. [21j tells us that x ~ — O.OGlGeV^ -3 which gives a negligible contribution 
to the vector meson mass close to chiral limit. 



CONCLUSION 

Based on the symmetries of the mixed (as well as unmixed) lattice QCD Lagrangian 
with Wilson and Ginsparg- Wilson fermions, heavy vector meson, partially quenched chiral 
perturbation theory was derived with explicit breaking of continuum symmetry up to order 
0(a 2 ). Using this formalism for the two flavor case, we extracted chiral and continuum 
expansion formulas with model dependent inputs with further comparison to the recently 
published numerical simulation which incorporated data over a wide range of meson masses. 

We perform our calculations in the isospin symmetry limit when m x = m y and m u = ma- 
in our calculations of vector meson masses we include the phenomenological contribution 
from the p — > tttt decay which can't be consistently extracted from heavy meson ChPT 
(because this is strongly relativistic process). We also didn't include the finite volume 
effects from the lattice. 



The comparison with the extrapolation expression successfully used in Ref. |21| leads us 
to the conclusion that in our approximation chiral and continuum symmetries are mostly 
broken by the operators in L mass and L a 2 and it is consistent to disregard the contributions 
from the Lagrangians: L a , L m 2 and L amq . This also gives us a simplified chiral-continuum 
extrapolation formula for further simulations using Wilson sea and Ginsparg- Wilson valence 
quarks. 

This study is important because it allows us to check the validity of lattice QCD simula- 
tions, QCD and ChPT as well, by performing systematic chiral and continuum expansions 
of lattice data to the real world with further comparison to the experimental data. 



a 4 + Xi m u = can be achieved if an = \2 = because m u is a variable. 
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